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CN , Abstract 

^ ! The breathing-mode giant monopole resonance is studied within the framework of the rel- 

ativistic mean-field (RMF) theory. Using a broad range of parameter sets, an analysis of 

constrained incompressibility and excitation energy of isoscalar monopole states in finite 

^ I nuclei is performed. It is shown that the non-linear scalar self-interaction considerably. It 

is observed that dynamical surface properties respond differently in the RMF theory than 

in the Skyrme approach. A comparison is made with the incompressibility derived from 

^ ■ the semi-infinite nuclear matter and with constrained nonrelativistic Skyrme Hartree-Fock 

S ' calculations. 
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The breathing-mode giant monopole resonance (GMR) has been a matter of contention in 
the recent past (Sharma 1990, Blaizot 1990). The energy of the GMR has been considered 
to be a source of information on the nuclear matter compressibihty. Theoretically, the 
incompressibility (compression modulus) is understood to have been obtained (Blaizot 1980) 
at about 210 MeV using the density-dependent Skyrme interactions. The deductions base 
themselves upon interpolation between various Skyrme forces for the GMR energies obtained 
from HF+RPA calculations. This approach has, however, not been able to reproduce the 
empirical GMR energies of medium- heavy nuclei. This is due to the reason that within 
the Skyrme approach (Treiner et al 1981) the magnitude of the surface incompressibility is 
about the same as the bulk incompressibility. 

Experimentally, the GMR has been observed in many nuclei all over the periodic table. 
However, precision data on the GMR have been obtained (Sharma et al 1988) only on a 
few nuclei including Sn and Sm isotopes. Analysis of the data (Sharma et al. 1989) using 
a leptodermous expansion led to the incompressibility of the nuclear matter at 300 ± 25 
MeV. This analysis assumed the third derivative of the EOS, which occurs in the Coulomb 
term in the expansion of the finite nuclear incompressibility (see Sharma 1990 for details), 
based upon the systematics from the Skyrme forces. In the Skyrme approach (Treiner et 
al 1981), the third derivative is related to the second derivative (incompressibility) of the 
equation of state (EOS). However, circumventing this relationship based upon the Skyrme 
theory, the reanalysis (Sharma 1991) of the data result in the incompressibility of the nuclear 
matter at slightly higher than 300 MeV and with increased error bars at about ±40 MeV. 
The data also allows to determine (Sharma 1991) the third derivative of the EOS. The 
results indicate that the magnitude of the surface incompressibility is much larger than the 
bulk incompressibility, in contrast with the relationship in the Skyrme ansatz. A detailed 
analysis of the GMR leading to the incompressibility of nuclear matter is, however, still 
under investigation. 

The relativistic mean-field (RMF) theory has achieved a considerable success (Gambhir, 
Ring and Thimet 1990, Reinhard 1989, Sharma, Nagarajan and Ring 1993) in describing the 
ground-state properties of nuclei at and far away from the stability line. On dynamical as- 
pects, however, only a few investigations have been made. The breathing-mode energies and 
incompressibilities were obtained (Maruyama and Suzuki 1989, Boersma et al 1991) within 
the RMF theory using the linear Walecka model. It is known that the surface properties 
of nuclei can not be described within the linear model. The correct description of the the 
surface properties and thus the overall properties of nuclei requires the inclusion of the non- 
linear scalar self-interaction (Boguta and Bodmer 1977). In the RMF theory the dependence 
of GMR energies on the surface properties of nuclei and on the incompressibility of nuclear 
matter is not yet known. In the non-relativistic Skyrme approach, on the other hand, the 



GMR energies are connected to the incompressibility of the bulk and the surface in a simple 
way. In this letter, we investigate the dependence of the breathing-mode energies in finite 
nuclei on the incompressibility of nuclear matter in the RMF theory and study the influence 
of the non-linear scalar self-interaction and surface properties on the GMR energies. 

We start from a relativistic Lagrangian (Serot and Walecka 1986) which treats nucleons 
as Dirac spinors ip interacting by exchange of several mesons: scalar a-mesons that produce 
a strong attraction, isoscalar vector u;-mesons that cause a strong repulsion and isovec- 
tor p-mesons required to describe the isospin asymmetry. Photons provide the necessary 
electromagnetic interaction. The model Lagrangian density is: 
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where U{cr) is the non-linear scalar self-interaction (Boguta and Bodmer 1977) with the 
cubic and quartic terms, required for an adequate description of the surface properties: 

U{a) = ^m^a^ + ^g^a^ + ^gsa'. (2) 

M, rricr, m^ and rup are the nucleon, the a-, the u;-, and the p-meson masses, respectively, 
and (7ct, g^i^ gp and e^/47r=l/137 are the coupling constants for the a-, the u;-, the p-mesons 
and for the photon. The field tensors for the vector mesons are Vt^'^ = d^uji, — d'^Up, R^'^ = 
d^Pu — d^pp — gp{p^ X Pj,) and for the electromagnetic field F'"'^ = d'^A^, — d" Ap. The 
variational principle leads to the stationary Dirac equation with the single-particle energies 
as eigenvalues, 

hoi^iix) = eiiJi{x), (3) 

where 

ho = -ia. V +P{M + gMr)) + guUJ%r) + gpTsp^r) + e^^A%r). (4) 

Solving these equations self-consistently one obtains the nuclear ground state in terms of 
the solution ipi. 

In order to obtain the isoscalar monopole states in nuclei we perform constrained RMF 
calculations solving the Dirac equation 

(ho - Xr'^)ipi{x,X) =eiilJi{x,X), (5) 



for different values of the Lagrange multiplier A which keeps the nuclear rms radius fixed at 
its particular value 

(If 1^/2 

R=[jjr'p,{r)d\j , (6) 

where 

A 

Px{r) = pv{r, A) = ^ iljj{x, A)^i(x, A) (7) 

is the local baryon density determined by the solution ■^i(x, A). The total energy of the 
constrained system 

ErmfW = EjiMF['4'i{x, A)], (8) 

is a function of A (or the rms radius R) which has a minimum, the ground state energy 
^RMF — Ermf{0), at a = corresponding to the ground-state rms radius R^. 

This behaviour of the constrained energy (8) as a function of A (or R) allows us to 
examine the isoscalar monopole motion of a nucleus as harmonic (breathing) vibrations 
changing the rms radius R around its ground-state value Rq- Considering s = {R/Rq — 1) 
as a dynamical collective variable and expanding (8) around the ground-state point s=0 (or 
A = 0) we obtain in the harmonic approximation, 

ErmfW = El^p + ]^AKc{A)s\ (9) 

where 

is the constrained incompressibility (Jennings and Jackson 1980) of the finite nucleus. The 
second term in eq.(9) represents the restoring force of the monopole vibration. In order to 
obtain its associated inertial parameter we apply the method (Maruyama and Suzuki 1989) 
of making a local Lorentz boost on the constrained spinors ipi{x, A): 

i,.{x,t) = ^S{v)A{x,X), (11) 
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where S{v) is the Hermitian local Lorentz boost operator 



^(t7) = cosh[^] +a.^sinh[^], (12) 



with = tanh ^{\v\) and 7 = \/l — tp. The velocity field v = v{r,t) is then obtained by 
the continuity equation 

- f^^lix, A)^,(x, A) + V-E4(a:, A)a^,(x, A) = 0, (13) 

which, using eqs.(7) and (11), transforms into the form 

. dp{r,X) 



ds 



+ V-(p(r,A)t;(r,t)) = 0. (14) 



Up to first order in s the velocity field is determined by eq.(14) into the form v = ■~su{r)f/\f\ 
where 

The inertial parameter of the monopole vibration is then obtained as (Maruyama and Suzuki 
1989) 

BreiiA) = A-' J u\r) HRMFir) d\ (16) 

where u{r) is the velocity function (15) at A = and 'Hrmf{'^) is the Hamiltonian density. 
We have obtained Kc{A) from eq.(lO), Brei{A) from eq.(16) and the frequency of the 
isoscalar monopole vibration as 
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Kci^ (17) 

\Brel{A) ^ '' 



for a number of spherical nuclei. Various parameter sets such as NLl (Reinhard 1989), 
NL-SH (Sharma et al 1993), NL2 (Lee et al 1986), HS (Horowitz and Serot 1981) and LI 
(Lee et al 1986) with values of the nuclear matter incompressibility K]^m =211.7, 354.9, 
399.2, 545 and 626.3 MeV, respectively, have been employed in the calculations. The last 
two sets, HS and LI, correspond to the linear model without the self-coupling of the a- 
field. In addition, the set LI excludes the contribution from the p- field. The sets NLl, 
NL-SH and NL2 correspond to the non-linear model, each having a different nuclear matter 
incompressibility and surface properties. Moreover, NL2 is the only set having a positive 
coupling constant (73, Eq (2). Whereas the sets NLl and NL2 reproduce the ground-state 
properties of nuclei only close to the stability line due to the very large asymmetry energy, 
the set NL-SH describes also nuclei very far away from the stability line (Sharma et al 1994). 



We show the results of our constrained RMF calculations in Figs. 1 and 2, where the 
collective mass Brei{A), the constrained incompressibility Kc{A) and the associated excita- 
tion energies ujc for a few nuclei have been displayed. First, we consider Kc{A) in Fig. 1 
(a). The incompressibility of nuclei shows a strong dependence on the nuclear matter incom- 
pressibility K]\fM, with a few exception for light nuclei. For the linear force HS, Kc shows a 
slight dip from the increasing trend for ^°^Pb and ^°Zr, whereas for light nuclei '^'^Ca and ^^O 
the HS values are even smaller than the NL2 values. The dependence of the imcompress- 
ibility K{A) of finite nuclei on Knm in non-relativistic Skyrme calculations is different: it 
increases monotonically with K]^m (Blaizot 1980). This difference can be understood from 
the difference in the behaviour of the surface incompressibility in the two methods. In the 
Skyrme approach, the surface incompressibility has been shown to be Ks ~ —K^^m for the 
standard Skyrme forces (Treiner et al 1981). However, this does not seem to be the case for 
the RMF theory as shown by the dip at the HS values. Thus, the surface incompressibility is 
not necessarily a simple function of the nuclear matter incompressibilty in the RMF theory. 

Since the GMR energy depends strongly upon the inertial mass parameter, we next 
examine the mass parameter B.f.f.i{A) in Fig. 1 (b). The collective mass obtained in Eq. 
(16) from the velocity field shows an interesting dependence on the mass of the nucleus 
considered. For the heavy nucleus ^"^^Pb, the collective mass is about 0.9 MeV~^ and for 
^°Zr it is 0.5 MeV^^. For lighter nuclei it decreases to about 0.3 MeV^^, as shown in Table 
1. Except for LI, where there is no p- meson coupling, it does not, however, depend so much 
upon the parameter set used and therefore shows only little sensitivity to K]^m- This implies 
that the Knm dependence of the energy is then predominantly due to the Knm dependence 
of the incompressibility Kc{A). 

We now consider two realistic parameter sets NLl and NL-SH and compare their results 
in Table 1. The constrained incompressibilities Kc for NL-SH are higher than for NLl as 
discussed above in Fig. 1(a). It is interesting to compare the relativistic collective mass 
Brei{A) with the expression usually employed in the non-relativistic sum- rule approach, 

Bsr{A) = MRl (18) 

where the ground-state rms radius Rq is used. The ratio Bsr{A)/Brei{A) is shown in Table 
1 where the energy uc, eq.(17) and its approximation 
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have also been compared for the sets NLl and NL-SH. The observation that the collective 
mass Brei{A) decreases significantly (down to 50 %) for light nuclei, provides a hint for the 
influence of the surface in the value Brei{A). For heavy nuclei the approximate value Bsr{A) 
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is rather close to Brei{A). Consequently, the monopole energies uc and ujsr differ by about 
1 MeV for heavy nuclei and up to about 5 MeV for the lighter ones. 

The energy Ugr corresponds to the monopole excitation energy Ei usually obtained from 
nonrelativistic Skyrme forces within the sum rule approach (Bohigas et al 1979) and in 
nonrelativistic constrained Hartree-Fock calculations. In Table 1 energies uc and oOsr have 
been compared with such nonrelativistic constrained HF results Usky obtained from Eq. 
(19) with the Skyrme-type forces SkM and SIII having about the same nuclear matter 
incompressibility Ki>4m as the sets NLl and NL-SH, respectively. From Table 1 it can be 
seen that the Skyrme results (Treiner et al 1981) ujsky are actually close to the energies uJsr 
and differ significantly from the values oiuJc- This difference in the RMF constrained energy 
ujc from the Skyrme constrained energy uJsky is small for heavy nuclei. It, however, increases 
for lighter nuclei, where RMF shows smaller values. The difference is particularly significant 
for the higher incompressibility forces NL-SH and SIII. It arises naturally from the lower 
values of the collective mass Bj.ei{A) for light nuclei in the RMF theory as discussed above. 
Thus, the masses B^ei and Bsr have different dependences on surface. 

Fig. 2 shows the constrained breathing-mode energy ooc for different nuclei. For heavier 
nuclei ^°^Pb and ^°Zr, uc shows a behaviour similar to that shown by Kq in Fig. 1 (a), 
as the collective mass does not show much change with incompressibility (Fig. l.b). The 
stagnation in energy at NL2 and HS is reminiscent of the effect of surface compression as in 
Fig. 1 (a), thus reflecting the role played by the surface in the RMF theory. In the lighter 
nuclei this effect is even more transparent as uoc shows a dramatic decrease for HS. Thus, for 
^°Ca and ^^O, ooc is not related to K^m in a simple way because of the combined effect of 
Kc and the collective mass -Brei) where there is a decrease in ooc with K^m- Employing only 
schematic parameter sets, it has been shown (Sharma, Nagarajan and Ring 1994) earlier that 
within the relativistic quantum hadrodynamics the surface compression responds differently 
than in the Skyrme ansatz. 

A comparison of the empirical values of the GMR energies with ujc is worthwhile. The 
empirical values for ^°^Pb and ^"^Zr are 13.9±0.3 and 16.4±0.4 MeV, respectively. For the 
lighter nuclei the values are very uncertain. Systematics of the values for ^'^^Pb show that 
the empirical value is encompassed by the constrained calculations curve from K]^m = 280 
- 320 MeV. The ^°Zr value is, however, overestimated by the corresponding curve by about 
1.5-2 MeV. 

We have also carried out calculations of the incompressibility for flnite nuclei by scaling 
the density using semi-inflnite nuclear matter. The incompressibility Ka for a flnite nucleus 
can be written as: 

Ka = K^m + KsA-''' + kJ^^]\1;-{i-27^]. (20) 

V A / 5 ro \ Knm 7 



where major contribution to Ka arises from the volume {K^m) and the surface terms {Ks)- 
The surface incompressibihty Ks has been obtained by calculating the second derivative of 
the surface tension {a) with respect to the density for each change in the scaled density 
(Blaizot 1980) 

Ks = 47rro^{22a(po) + 9p^^"(po) + ^^^P^e'^(po)}, (21) 

The details of the procedure to perform scaling of the density have been discussed by Stocker 
and Sharma (1991). The last two terms in eq. (20) contribute very little. We add these terms 
for completeness, however. The asymmetry coefficient K^, has been taken at —300 MeV 
from the empirical determination (Sharma et al 1989) and is a reasonable value. The third 
derivative of the EOS, Pqc'", for each force is obtained from the nuclear matter calculations. 

We show the 'scaling' incompressibihty Ka obtained from Eq. (20) for various RMF 
forces in Fig. 3. The general trend of the scaling incompressibihty with K^m is about the 
same as in Fig. 1(a), showing a dip for the force HS. In general, the Ka values are about 20% 
larger than the constrained values Kc- This is consistent with the known relationship of 
the two types of the incompressibilities for nuclear matter, whereby it was shown (Jennings 
and Jackson 1980) that Kscai{NM) = j^Kconstr{NM). It is seen from Fig. 3 that from 
NLl to NL-SH, Ka values increase. However, the HS values show a clear deviation from 
the increasing trend, and particularly for the light nuclei where surface compressibility and 
surface properties play increasingly significant role. Thus, Fig. 3 brings out again the 
importance of the role of the surface in the dynamical calculations within the RMF theory. 
The surface incompressibihty Ks for the forces NLl and NL-SH are —333.1 and —610.1 
MeV respectively. Here the ratio of the surface to the bulk incompressibihty is obtained as 
1.58 and 1.72 respectively. These values differ considerably from the ratio of 1 (Treiner et 
al. 1981, Sharma et al. 1989) in the Skyrme ansatz. The values of the bulk and surface 
incompressibilities for NL-SH are closer to the empirical values from Sharma et al (1988). 

In conclusion, it has been shown that within the RMF theory, the incompressibihty of a 
finite nucleus, K{A), does not depend on the nuclear matter incompressibihty in a simple 
way and that it is influenced significantly by the properties of the surface subjected to 
compression. This admits the possibility to obtain the surface incompressibihty of nuclei 
with a magnitude different from that of the bulk incompressibihty. This behaviour of the 
surface in the RMF theory is different than in the Skyrme ansatz. The collective mass for 
the breathing mode vibration in the RMF theory also shows a different behavior from light 
to heavy nuclei, thus affecting the frequency of the isoscalar breathing mode. 
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TABLES 
?ABLE I. The constrained incompressibility Kc{A) in MeV, the mass parameter Brei{A) in 

'^^, the ratio Bsr{A)/Brei{A) and the associated monopole frequencies loq and ujsr, (both in 

') calculated with the sets NLl and NL-SH. Comparison is made with the constrained Skyrme 

ts uJsky obtained within the sum-rule approach (Bohigas et al 1979) with the Skyrme forces 

and SIII. 

NLl: Knm = 211.7 MeV NL-SH: Knm = 354.95 MeV 

SkM: Knm = 216.7 MeV SIII: Knm = 356.00 MeV 



ei Kc{A) B{A) Bgr/Brel UJC ^sr ^Sky Kc{A) B{A) Bsr/Brel UJC ^sr ^Sky 

) 74.0 0.2282 0.742 18.0 20.9 22.4 105.6 0.2911 0.544 19.0 25.8 26.6 

a 102.0 0.3086 0.894 18.2 19.2 20.2 153.4 0.3578 0.750 20.7 23.9 24.7 

r 117.1 0.4776 0.924 16.2 16.3 17.0 149.0 0.4677 0.930 20.4 18.5 21.2 

'6 116.0 0.7825 0.991 12.2 12.2 12.9 197.7 0.8084 0.939 15.6 16.1 16.2 



FIGURE CAPTIONS 

xre 1. (a) The constrained incompressibility Kq in Eq. (10) for a few nuclei 
ined using various RMF parameter sets, (b) The collective mass B^ei for the 
thing-mode monopole vibrations obtained from Eq. (16) within the RMF theory. 

ire 2. The frequency ojc of the monopole mode obtained using eq. (17). 

are 3. The incompressibility Ka (Eq. 20) obtained from 'scaling' of the nuclear 
lity in the semi-infinite nuclear matter using the Thomas-Fermi approximation. 
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